Ground-State Entanglement in Interacting Bosonic Graphs 
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We consider a collection of bosonic modes corresponding to the vertices, of a graph r. Quantum tunneling can 
occur only along the edges of F and a local self-interaction term is present. Quantum entanglement of one vertex 
with respect the rest of the graph is analyzed in the ground- state of the system as a function of the tunneling 
amplitude r. The topology of F plays a major role in determining the tunneling amplitude r* which leads to 
the maximum ground- state entanglement. Whereas in most of the cases one finds the intuitively expected result 
r* = oo we show that it there exists a family of graphs for which the optimal value of r is pushed down to a 
finite value. We also show that, for complete garphs, our bi-partite entanglement provides useful insights in the 
analysis of the cross-over between insulating and superfluid ground states 

PACS numbers: 03.67, 03.67.L 



Introduction- Entanglement measures quantify the strength 
of purely quantum correlations between subsystems of a com- 
pound quantum system 1 1 ] . In the last few years efforts in the 
new field of quantum information science unveiled how such 
correlations can be exploited as a genuine resource for carry- 
ing out computational and communication tasks beyond the 
reach of any classicaly operating device |2]. More recently 
several studies pointed out that the notion of quantum entan- 
glement can be a useful conceptual tool to investigate the com- 
plex properties of quantum many-body systems; in particular 
the role of entanglement has been analyzed in spin systems 
undergoing quantum i.e., ground- state phase transitions 1 3] 

In this paper we shall study a related problem: the entan- 
glement behaviour in the ground-state of a system made of a 
finite number of bosonic modes bi-linearly coupled each-other 
and with a repulsive self-interaction. Such a system can by 
described by a graph 01 whose vertices are associated with 
the bosonic modes themselves and whose edges correspond 
to the bilinear couplings i.e., tunneling. The corresponding 
quantum Hamiltonian is a Bose-Hubbard one 01 that in the 
general case represents a very difficult many-body problem. 
In order to effectively tackle this problem we will mostly fo- 
cus on graphs with a small number of vertices i.e., four. In 
spite of this simplification our analysis reveals a variety of 
features that are expected to be of general validity. 

In Ref 1 6] it has been shown that, in absence of of non-linear 
self-interaction ground- state entanglement of one vertex with 
respect the rest of the lattice contains information about the 
graph topology. Moreover the graph topology has a subtle in- 
terplay with self-interaction in affecting the entangling power 
01 of tunneling coupling in small graphs 01 . 

In the following we will adhere to the approach to quantum 
entanglement i n sy stems of indistinguishable particles dis- 
cussed in Refs |l d|llLll3.ll3ll . A complementary one is pur- 
sued in Refs. !14ll5Lll6 1 and more recently in 1 17]. In our ap- 
proach the subsystems are provided by bosonic modes and not 
by the particles; in particular even the notion of single parti- 
cle entanglement makes sense. This mode-entanglement con- 
cept is an inherently second-quantized one: the tensor prod- 
uct structure of the state-space necessary in order to define 
entanglement 1 18] is provided by identification of the bosonic 



(fermionic) Fock space with a set of linear oscillators (qubits) 
associated to single-particle state vectors. 

As it will be illustrated in the last part of the paper the bi- 
partite mode entanglement we study is closely connected to 
local i.e., on-site, particle number variance. Quite recently 
this quantity has been showed to be a truy physical resource 
to overcome limitations imposed by mass superselection con- 
straints by teleportation 1 19]. Another quantum-information 
theoretic motivation to our work is to study how interaction 
and graph topology affect the capability of entanglement gen- 
eration i.e., the entangling power 1 7],| 8], by adiabatically turn- 
ing on the tunneling parameter at zero-temperature 

Preliminaries- By an interacting bosonic graph we mean 
a collection of bosonic modes associated to the vertices of a 
graph T := (V,E) [4] whose dynamics is governed by the 
following Bose-Hubbard hamiltonian Q] 



H = 



• (b\b j +h.c.)+e y E t 
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where rii := b\bi is the occupation number operator of the 
i-th vertex of F. Eq. {U can describe an inomogenehous spa- 
tial structure in which bosonic particles e.g., ultra-cold atoms, 
reside over a collection of sites -the set V of vertices of T- 
and can tunnel among different locations, with amplitude r, 
across the edges E of the graph. The non-linear part in Eq 
{T}, wheighted by the interaction parameter e, accounts for 
the self-interaction of the modes when more than one particle 
is present in the same site. 

It is important to stress that the Hamiltonian |l} may de- 
scribe a variety of quantum physical systems: photonic modes 
coupled by beam-splitters and with Kerr non-linearities, ar- 
rays of Josephson junctions, ultra-cold atoms in some (inomo- 
genehous) optical lattice B2QI1 . For the sake of concreteness 
we will mostly use a language in which the bosonic modes 
are thought to be spatially localized e.g., optical-lattice sites. 
These kind of systems have been already considered in the 
quantum information literature in Refs. 12(1 EH El El E3. 

13 His. 

In ref. the ground- state entanglement properties have 
been analyzed for pure tunneling. The main purpose of this 
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work is to extend those results to r < oo. In this regime the 
dynamics described by {T} is a complex one due to the pres- 
ence of two competing effects. The tunneling term in JO is re- 
sponsible for delocalizing the particles over the graph vertices, 
whereas the non-linear coupling term tends to localize them. 
This interplay is dramatically displayed by the occurence of 
the superfluid-insulator transition predicted by the Hamilto- 
nian ([0 over lattices with a large number of sites |5] for a 
critical value of the coupling strength e. One is then naturally 
led to conjecture that the higher the tunneling amplitude r the 
greater a single mode gets entangled with the rest of the graph. 
We will see that that turns out to be the case for many of the 
possible graph topologies but that it also exist topologies for 
which an increasing of r results in an decreasing of the mode 
entanglement of one vertex. 

Mode-entanglement and self-interaction- In this section we 
give a qualitative description of the results of the computations 
for the different systems under study. Our goal is to charac- 
terize, for the reference mode 0, the behaviour of the mode 
entanglement of the ground state \S&gs(t)) of Hamiltonians 
of the kind |l} when the ratio between the hopping param- 
eter and the self-interaction parameter r/e varies from zero 
to a value that is much greater then one. In order to do so 
we fix e = 1, that is we measure r in e units and let r vary 
in [0, Tm ax ] (r max > 0) with step At. For each value of r 
we compute E(r) = Tr[po(r) log po(r)]/log(N + 1), where 
po(r) = Tr Vxo (|^M)(^(r)|) i.e., the reduced density ma- 
trix of the mode HI, the logarithms are taken in base two 
and the factor 1/ log(7V + 1) gives the normalization of the 
mode entanglement to its maximum possible value. In our 
first simulations we have considered rooted graphs Tj with 
N = L = 4, see figure (1); their Hamiltonian is given by 
where the tunneling is allowed only between the sites linked 
by the edges of the relative graph. 

We first focus on the interval in which r < 1. An interst- 
ing feature that can be highlighted in this regime, in which the 
tunneling part of the Hamiltonian can be considered a pertur- 
bation of the self-interaction part, is the ordering of the curves 
with respect to the graph topology. This feature is particu- 
larly clear in fig. (4) where the graphs Tj, j = 11, 12, 13 can 
be obtained by adding respectively 1 , 2 and 3 links to the sub- 
graph of Tio. One can see that the greater the connectivity of 
the sub-graph Tj — {0} the greater the entanglement, that is 
Er 13 > Er 12 > Er 1± > Er 10 . The same ordering appears for 
the other two sets of graphs. One finds Er 5 > Er 3 > Er 4 , 
where the sub-graphs of T% and T4 have both two links; while, 
see figure (3), Er 9 > Er 7 > Er 8 > Er 6 , and in this case the 
sub-graphs of T7 and have the same number of links i.e, 
two. The ordering of the graphs according their ground- state 
entanglement for small tunneling coupling r can be related to 
the spectrum of the one-particle tunneling Hamiltoanian i.e, 
the spectrum of the adiacency matrix of the graph. Indeed by 
diagonalizing the adiacency matrix of all the L = 4 graphs we 
found that the greater (in modulus) the minimum eigenvalue 
So the greater the associated entanglement. This seems to be 
a natural and general result. Indeed \eo\ gives the strength of 
the tunneling rate in the ground- state with no self-interaction, 
once r is turned on, one expectes the full ground-state to be 



mostly a coherent mixing of the the ground state with r = 
i.e., the state with one particle per vertex, and the pure tun- 
neling ground-state. The greater so the greater the wheight of 
this latter. 

Noticeably in the regime r >> 1 the ordering of the curves 
is inverted. This feature is again very clear for the set of 
graphs Tj, j = 10, 11, 12, 13; in fact, as we can see in fig. (4), 
Er 10 > Er 1± > Er 12 > Er 13 i.e., the greater the connec- 
tivity the lower the mode entanglement of the mode 0. This 
behaviour starts to be evident when r > 1 and is mantained 
even in the asymptotic regime (r ^> 1) where now it is the 
self-interacion part of the Hamiltonian that plays the role of 
the perturbation. The latter feature can be seen by looking 
at the inset of figure (4) where it is displayed the behaviour 
of E(r) over the full interval [0,r ma J, r ma;E = 20; when 
r ^> 1 the ordering of the curves remains the same described 
for r > 1. For the set of graphs described in figure (3) we 
have the same kind of behaviour: the graph with the lower 
connectivity ^(one link in the sub-graph) presents the high- 
est value of E(t), while Tg, the graph with the higher connec- 
tivity (three links in the sub-graph) displays the lowest val- 
ues of E(r). In this case of T? and T$ the ordering in the 
asymptotic regime remains the same displayed for r <C 1 i.e., 
Er 7 > Er 8 . 

We treat separately the results (see figure 2) for the graphs 
of the set a) because in this case an interesting behaviour 
comes into play. Here the ordering of the curves for r ^> 1 is 
of the same kind of the one seen for the other set of graphs: 
T$, the graph with the higher connectivity (three links in the 
sub-graph) has the lowest values of E(r). But the very in- 
teresting feature is that, whereas for all the other graphs we 
have seen so far E(r) grows monotonically as a function of 
r, in the case of T4 and T$ this is no longer true. In fact 
what happens is that when r starts to be different from zero 
E(r) grows in both cases but then for higher values of r it be- 
comes a decreasing function and reaches a stationary value for 
r ^> 1. This behaviour is more evident for Ts, for which E(r) 
presents a maximum for r = 1.14, but it is also characteristic 
of T4, for which E(r) presents a maximum for r = 3.22. 
This peculiar behaviour can be somehow understood analyt- 
ically in the following way. Let us consider for example the 
system represented by T$; the latter belongs to a class of sys- 
tems whose Hamiltonian, in the pure tunneling regime (e = 
0), can be written as H = — t(6q6i +h.c.) — r Yli^j=i b\bj. It 
corresponds to the topology T in which the vertex is con- 
nected only with the vertex 1 whereas the subgraph with ver- 
tices {1, . . . , L} is a complete one. Now we provide a simple 
argument showing that, in the large tunneling amplitude limit 
and for a large number of vertices, the ground- state entangle- 
ment of the vertex for the systems described by the above 
Hamiltonian is vanishing. 

By introducing the Fourier operators bk '•= 
l/VL^TE^=i 1 exp(^)6 J (k = 0,...,L - 2) 
the Hamiltonian associated to T 5 (with e = 0) can 
be rewritten as H = —r(L — I) (Ho + Hi), where 
H = 1%, Hi = (L — l)- 1 [(6jfei + hx.) - JV L _i],where 
Nl-i := Ylf=i bjbj. The second term in the equation above 
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FIG. 1: Inequivalent rooted graphs for iV = L = 4; the root vertex 
is directly linked with 1 (a), 2 (b) and 3 (c) vertices. 




FIG. 2: Plot of E(t) for the systems T jf j = 3, 4, 5. 




FIG. 3: Plot of E(r) for the systems T h j = 6, 7, 8, 9. 
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FIG. 4: Plot of E(t) for the systems T jf j = 10, 11, 12, 13. 



can be regarded, for L oo, as a small perturbation of the 
first with coupling constant ~ It follows that the ground 
state of Ho with N = L particles, for L — oo is given by a 
condensate over the mode bo, i-e., \^gs)l=oo ~ (^o) L |^)- 
In this latter state the mode clearly factors out with zero 
occupation number, therefore we have vanishing mode 
entanglement. A simple first-order pertubation evaluation 

gives \^ gs )l ~ \^gs)l=oo - L-^Jft^-^O) which 
shows that, for large enough L, the entanglement of the zero 
mode is a monotonically decreasing function of L. This large 
L behaviour is anticipated by the system T 5 : here L is finite 
and small, therefore E(r) does not go to zero but it decreases 
and it reaches a finite non-zero value in the asymptotic regime 
i.e., r ^> 1. The result presented above is robust against the 
turning on of a small self-interaction coupling e ^ Indeed 
even such a terms would be order L _1 with respect Ho ll3lll 

This discussion indicates that the size of lowest single parti- 
cle eigenvalue of T — {0} plays a major role in determining the 
entanglement properties we analyze in this paper. This quan- 
tity in turn is well-known to have a clear topological meaning, 
for example for a regular graph with order r vertices one has 
So — —r [4]. Roughly speaking, the greater the connectivity 
of T — {0} the greater \eo \ • 

Insultator-Super fluid cross-over- In this section we show 
that the kind of bi-partite entanglement we analyzed in this 
paper provides useful insights on the itinerant vs localized 
character of the particles in the ground state. In particu- 
lar E(t) can be related to the local particle number vari- 
ance that it is a standard tool to study the insulator- superfluid 
transition [30]. We give a description of the simulation re- 
sults of for the case of complete graphs. To start with one 
can consider the simplest possible case, given by L = 2 
i.e., the bosonic dimer 112 611. It is elementary to see that 
the ground state of 0} with e = 1, given then given by 
|*gs(t)> := cos(0/2)|ll> + sin(0/2)(|O2> + |20»/>/2, 
where 6 = — tan -1 (2r/e) from which it follows immedi- 
ately that the ground- state entanglement is given by E(r) — 
-cos 2 ((9/2) In cos 2 ((9/2) - sin 2 (<9/2) ln(sin 2 (<9/2)/2). It is 
imeediate to see that E(r) is a monotonic increasing function 
of r; maximal entanglement is then achieved in the pure tun- 
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FIG. 5: Plot of : a) E{t), b) derivative of i2(r), c) derivative of 
(nf ) — (rii) 2 for the complete graphs with j = 3, 4, 5, 6, 7 

neling regime. 

In order to measure the itinerant character of the particles 
over the graph is useful to analyze the ground- state variance 
of the local occupation numbers i.e., (nf ) — (n^) 2 , (i E V). 
This quantity plays the role of a sort of order parameter in 
the insulator- superfluid transition: small (large) values of it 
are associated to an insulator (superfluid). In t he dimer case 
one gets the variance sin 2 (<9) = 1/2(1 - 1/Vl + 4r 2 ) that 
it is again a monotonic increasing function of r; its derivative 
shows a peak for r = l/2y/2. The same kind of qualitative 
behaviour is displayed by E(r) and by its derivative. 

We considered complete graphs corresponding an increas- 
ing number of sites to i.e., N = L = 3,4,5,6,7, see fig- 
ure (5). In this case E(r), which is plotted in the graph- 
ics a), increses monotonically for all the interval [0, r max ], 
Tmax = 1. In the inset of the figure (5) it is plotted the vari- 
ance (n 2 ) — (n^) 2 , (i G V), that, in view of the symmetry of 
the systems, is the same for all the sites i. As for the mode 
entanglement it increases monotonically for all the interval 
[0, T max ]. We have therefore plotted the derivative of this 
quantity and we have compared it with the derivative of the 
mode entanglement, see graphics c) and b) in figure (5). Of 



course the insulator- superfluid phase transition occurs only in 
the thermodynamic limit (when A 7 ", L — > oo) so the peak in the 
derivative plot can only be regarded as a far precursor of this 
transition. In this case, the derivative exibits the behaviour 
expected by such a precursor, in fact, as TV = L increases 
it becomes more and more peaked; in presence of an actual 
phase transition it should display a singular behaviour in cor- 
respondance of a critical value r c of the tunneling parameter. 
These results clearly show that E{r) has the same qualitative 
behaviour of the variance the local particle variance [32] and 
therefore, at least for complete graphs, it represents on itself 
a useful quantity to analyze the cross-over between insulating 
and super-fluid phases. We observe that this result is remines- 
cent of an analog one for the metal-insulator transition in the 
Fermionic Hubbard nodel B33I1 . 

Conclusions - We have analyzed ground- state entangle- 
ment of interacting bosons for all the four vertices graphs. The 
ground-states for all possible graph topologies have been ob- 
tained by numerical diagonalization and the entanglement of 
one vertex with respect to all the others has been computed 
as a function of the tunneling amplitude. This analysis al- 
lows to order graphs, for a given tunneling amplitude r, in 
terms of the amount of this bi-partite ground- state entangle- 
ment E(t). Remarkably this order actually depends on both r 
and the graph topology: for small (large) r the higher (lower) 
the modulus of the minimum single-particle eigenvalue of the 
graph T (sub-graph T — {0}) the higher (lower) the mode en- 
tanglement. Moreover numerical results show un unexpected, 
non monotonic, behaviour of E(r) for some particular graph 
topologies. Such a phenomenon can be understood analyti- 
cally in the limit of no self-interaction. We finally showed, 
for complete graphs of different sizes, that E(r) contains use- 
ful information about the cross-over between the insulator and 
superfuild regime. The analysis of the role of topology and 
self-interaction on multi-partite entaglement, for example by 
considering all possible bi-partitions of the graph at once, is a 
more demanding task clearly deserving future investigations 
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